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We study the zero-bias magnetoresistance (MR) of an interacting quantum dot connected to 
two ferromagnetic leads and capacitively coupled to a gate voltage source Vg. We investigate the 
effects of the spin-activity of the contacts between the dot and the leads by introducing an effective 
exchange field in an Anderson model. This spin-activity makes easier negative MR effects, and can 
even lead to a giant MR effect with a sign tunable with Vg. Assuming a twofold orbital degeneracy, 
our approach allows to interpret in an interacting picture the MR(Vg) measured by S. Sahoo et al. 
[Nature Phys. 2, 99 (2005)] in single wall carbon nanotubes with ferromagnetic contacts. If this 
experiment is repeated on a larger V^g— range, we expect that the MR(V9) oscillations are not regular 
like in the presently available data, due to Coulomb interactions. 

PACS numbers: 73.23.-b, 75.75.+a, 85.75.-d 



I. INTRODUCTION 

The quantum mechanical spin degree of freedom is now 
widely exploited to control current transport in electronic 
devices. For instance, the readout of magnetic hard disks 
is based on the spin-valve effect, i.e. the tunability of 
a conductance through the relative orientation of some 
ferromagnetic polarizationsi. However, realizing spin in- 
jection in mesoscopic conductors would allow to imple- 
ment further functionalities, like e.g. a gate control of 
the spin valve effect^-'^. Importantly, electronic interac- 
tion effects can occur in mesoscopic structures, due to 
the electronic confinement. This raises the fundamental 
question of the interplay between spin-dependent trans- 
port and electronic interactions. 

Upon scattering on the interface between a ferromag- 
net (F) and a non-magnetic material, electrons with spin 
parallel or antiparallel to the magnetization of F can pick 
up different phase shifts, because they are affected by 
different scattering potentials. This Spin- Dependence of 
Interfacial Phase Shifts (SDIPS) can modify significantly 
the behavior of mesoscopic circuits. First, when a meso- 
scopic conductor is connected to several F leads with non 
collinear polarizations, the SDIPS produces an interfacial 
precession of spins which can modify current transport in 
the device^'S'Si^. Secondly, in collinear configurations, 
precession effects are not relevant, but the SDIPS can 
modify mesoscopic coherence effects. For instance, in su- 
perconducting/ferromagnetic hybrid circuits, the SDIPS 
introduces a phase shift between electron and holes corre- 
lated by Andreev reflection^. References ^3 and have 
identified signatures of this effect in the experiments of 
Refs. ^3 and 0, respectively. In principle, normal sys- 
tems in collinear configurations can also be affected by 
the SDIPS. Indeed, from Ref. Ti, the SDIPS should pro- 
duce a spin-splitting of the resonant states in a ballistic 
interactionless wire contacted with collinearly polarized 
ferromagnetic leads. However, this has not been con- 
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FIG. 1: Mesoscopic element M connected to ferromagnetic 
leads L and R. The magnetic polarizations pl and pr of 
leads L and R can be parallel (configuration P) or antiparallel 
(configuration AP). The element M is capacitively coupled to 
a gate voltage source Vg. 



firmed experimentally ye^i^. 

Recently, Ref. ^3 has reported current measurements 
in a single wall carbon nanotube (SWNT) connected to 
two ferromagnetic leads with collinear polarizations. The 
asymmetries observed in the magnetoresistance (MR) of 
the SWNT versus gate voltage are strikingly similar to 
those predicted by Ref. 14 for an interactionless wire 
subject to the SDIPS^*. However, the SWNT of Ref. ^ 
showed a quantum dot behavior with strong Coulomb 
Blockade effects, as demonstrated in a great number 
of experiments with non-magnetic leads (see e.gi^^). 
Therefore, one important question is how interaction 
effects modify the scheme proposed by Ref. 14. The 
problem of the effects of interactions on the transport 
properties of a central region connected to ferromag- 
netic contacts has already been considered in various 
regimes, like e.g. the Coulomb blockade regimaSiSiiSSiS^, 
the Kondo regim o^^i^^ , the Luttinger liquid regim o^i^^ 
and the marginal Fermi liquid regime^. This article de- 
velops an approach suitable for the limit of Ref. [13 and 
studies, for the first time, the effect of the SDIPS on a 
quantum dot. We consider a quantum dot coupled to 
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metallic leads through spin-active interfaces. We use an 
Anderson model to study the MR of the circuit above 
the Kondo temperature, but beyond the sequential tun- 
neling limit. The SDIPS is taken into account through 
an effective spin-splitting of the dot energy levels. This 
splitting makes easier negative MR effects. When it is 
strong enough, it can even lead to a giant MR with a 
sign oscillating with the dot gate voltage Vg , similarly to 
what has been found in the non-interacting case. In the 
non-interacting case, assuming that the properties of the 
contact are constant with energy and that the SDIPS is 
too weak to split the conductance peaks, one finds that 
the MR(Vg) pattern is similar for all conductance peaks. 
In contrast, the effect of the SDIPS depends on the occu- 
pation of the dot in the interacting case. This is in appar- 
ent contradiction with the data of Ref. ilTj because, in the 
V^— range presented in this Ref., the MR(V^) oscillations 
are regular. Using a two-orbitals model, which takes into 
account the K — K' orbital degeneracy commonly ob- 
served for SWNTs (see e.g. Refs. l2l2l29ii3dElli33) . one 
can solve this discrepancy. In this framework, we expect 
non-regular MR(V^) oscillations if the experiment is re- 
peated on a larger V^— range. 

This article is organized as follows: we start with sum- 
marizing the results found for the non-interacting case in 
section lll Al Then we introduce a model for the interact- 
ing case in section lll Bl Section lllll addresses the case of a 
one-orbital quantum dot circuit, and section Hvl the case 
of a two-degenerate-orbitals quantum dot circuit. Finally 
section concludes. 



the conductance G'^(Vg) of the circuit in configuration c 
presents Fabry-Perot like resonances for „ ~ 0, with 
= {27Td - - ^"AnC^TrN^)-' - eaVg - < a 
resonant energy, c? S Z, Ep^ the wire Fermi energy, 
the density of orbitals states at the Fermi level in the 
wire and a the spin direction opposite to a (we have 
used eaVg Ep). From this Eq., in configuration c, 
the SDIPS produces a spin-splitting 

gf^BhsDiPS = Ed,i - ^^,1 = Am '^^^ 

ie{L./?.} ^^^^ 

of the resonant energies. When the effective field h'^sDiPS 
is strong enough to produce a spin-splitting of the con- 
ductance peaks, the circuit can display a giant MR effect 
with a sign oscillating with Vg, due to the strong shift 
of the conductance peaks from the P to the AP con- 
figurations. In the opposite case, MR remains smaller, 
but the SDIPS can still be detected through character- 
istic asymmetries in the oscillations of MR versus Vg 
(see Fig. 2- right of Ref. IT^ . Importantly, assum- 
ing that Ti''^ and v?^''^ are constant with Vg, one has 
G^[Vg) = CiVg + [2/eaNlf]) . This implies that when 
hgpijpg is not strong enough to produce a spin-splitting 
of the conductance peaks, the MR(Vg) pattern is similar 
for all the peaks displayed by G^{Vg). 

B. Interacting case 



II. MODEL 

We consider a mesoscopic element M connected to fer- 
romagnetic leads L and R (Fig. The chemical poten- 
tial of M can be shifted by eaVg using the gate voltage 
Vg, with a the ratio between the gate capacitance and 
the total capacitance of M. The magnetic polarizations 
PL and pji of leads L and R can be parallel (configuration 
c = P) or antiparallel (configuration c = AP). 

A. Non-interacting case 

Before introducing the interacting model investigated 
in this article, it is useful to reconsider the results ob- 
tained by Ref. ^3 for the case in which M is a non- 
interacting single-channel ballistic wire of length L. In 
a scattering approach"^'^-, the conductance of the circuit 
depends on the transmission probability T^^''^ for elec- 
trons with spin cr £ {Ti i} through contact I G {L,R}, 
and on the reflection phase ip^''^ for electrons with spin 
a coming from the wire towards contact I. The in- 
dex c = P[AP] denotes the parallel [antiparallel] leads 
configuration. A spin dependence of tp'^''^ can occur 
due to the magnetic properties of the contact materi- 
als used to engineer lead I. Due to size quantization. 



We now assume the presence of strong Coulomb in- 
teractions inside M, such that we have a quantum 
dot connected to ferromagnetic leads. Such a sys- 
tem can be realized for instance by using granular 
films^, nanoparticles^^, carbon nanotubes^'^^''^®, or Cgo 
molecules^i. In the non-interacting case of section III Al 
we have considered that the spin-dependent confinement 
potential felt by electrons causes the SDIPS, which leads 
to the spin-splitting of the resonant states. In the in- 
teracting case, the scattering approach is not suitable 
anymore. However, the energy of the quasi-bound sin- 
gle particle states in quantum dot M can depend on spin 
due to the spin-dependent confinement potential. On this 
ground, we adopt the effective Anderson hamiltonian 

H — Hdot + Hleads + He (2) 

with 

d,a d,d' ,(7,(t' 

(d,a)it(d' 

Hleads — ^ ^ ^ka^j^fj^ka 
k,a 

= E [^dAa^^ka + (<L)*cLcd.) 
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Here, ^do- refers to the energy of the dot orbital state d 
for spin a, to the energy of lead state k for spin a and 
t^^ is an hoping matrix element (we assume that the spin 
a is preserved upon tunneling like in section FlI A|) . The 
index k runs over the electronic states of lead L and R. 
Coulomb interactions are taken into account through the 
term in U — e^/Cs, with nda — c\^Cda and Cs the total 
capacitance of the quantum dot M. By construction of 
the model (see above), for [/ = 0, each orbital level ^do- 
corresponds to a resonant level E'^^^ of section III Al with 
Cdi — CrfT = dt^Bh'^sDiPS- '"^'^ therefore regard the ef- 
fective Zeeman splitting hgjjjpg in model l(2Jl as a gener- 
alization of the SDIPS concept to the interacting caseS. 
The specificity of this effective field, with respect to an 
ordinary external field, is that it depends on the configu- 
ration c of the ferromagnetic electrodes. For instance, in 
the case of symmetric ferromagnetic contacts, symmetry 
considerations lead to h^jjips 7^ and h^DiPS ~ ^■ 

In the following, we calculate the zero-bias conduc- 
tance of the circuit usingS^ 



with 



h G" 



rL(fi^)rS.(M 



(3) 

Im[Gdcr(w) 



The above equation involves the retarded Green's func- 
tion GdA^) = j^^Gda{t)e"^'dt with Gda{t) = 

—iO{t)(JyCda{t)-,c\^{0)Yj. We also use the Fermi dis- 
tribution /(^) = (1 -|-exp[^])~^ and the tunnel transition 
rates T\^{0 = E 2^ ji^ l' '5(e = ik.) with / e {L,R}. 



E27r|tL| '5( 

fc 

Note that Gda, ^da- and F^^ depend on the configuration 
c £ {P, AP} considered but for simplicity we omit the 
index c in those quantities. We want to study current 
transport in the limit studied in Ref. UtI i.e. the width 
of conductance peaks displayed by the circuit is deter- 
mined not only by temperature but also by the tunnel 
rates (F^^. -I- F^^ ~ 2fcsT). This requires to go beyond 
the sequential tunneling description, i.e. to take into 
account high-order quantum tunneling processes. For 
this purpose, we will calculate Gda using the equation 
of motion (E.O.M.) technique^, which is valid for tem- 
peratures larger than the Kondo temperature Tk of the 
systemjii. 



III. SINGLE LEVEL QUANTUM DOT 

For simplicity, we first take into account a single orbital 
level d of the dot. We follow the lines of Ref. \M The 
E.O.M. technique leads to 



Gdaito) 



1 - {rida) 



yrida) 



(4) 



/ + 00 T 
— /(fi^)Im[G^(c^)] (5) 
-00 I" 

the average occupation of orbital d by electrons with spin 
a. We define 



S^' = J:"da~UJ:'d:,c^[h^-^da^U-J:''da-^ltd^r' : (6) 



U + Y. 



-\ (7) 



(8) 



and, for i G {1, 2, 3}, 



E 



Aii(6(T') 



fd 



h^-idc+ £.d'cr' - £.ka' + iO+ 



E 



kcr' 



^ - £.da - CdV - nil + £,ka' + iO+ 



(9) 



Here, one has /ii(C) = /(O, ^l2iO = 1 - fiO and 
fi^ — I (We anticipate on the next paragraphs by defining 
^dad'a' u & N and an arbitrary dot state d'a' ^ da, 
but only n — 1 and d'a' = da are needed for the present 
one-orbital case). We assume that the coupling to the 
leads is energy independent (broad band approximation), 
which gives e.g. 2°, = -z(Fi, F« )/2. The term 2°^, 
which is due to the tunneling of electrons with spin a, 
already occurred in the non-interacting case^. In the in- 
teracting case, Gdaii-^) also involves S^'"^,^, terms related 
to the tunneling of electrons with spin a. The average 
occupation (uda) can be calculated from Eqs. Q and 
(O as 



{nda) 



{nsa) (1 - {nsa)) + {nsa) (noa) 



1 - {{noa) - {nsa)){{nDa) - (nsa)) 
with, for j e {S,D}, 



■ Im 



+00 



hdhj 



^ - £,dcr - Sct(w) 



(10) 



(11) 



Figure [3 shows the conductance G^ in configuration 
c G {P, AP^ (panels a, c and e) and the magnetore- 
sistance MR = (G^ - G'^^)/(G^ + G^^) (panels b, d 
and f) calculated for different values of hf^sDiPS^ using 
F^^j^j — F;(l ±Pi) for I G {L, i?}. We have used parame- 
ters consistent with Ref. [TtI i.e. U/ksT = 30, and tunnel 
rates Yl(r) leading to the proper width and height for 
the conductance peaks. We have also used relatively low 
values for Pl(r) because usual ferromagnetic contact ma- 
terials are not fully polarized"*"^. The conductance peak 
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FIG. 3: Average occupations {nd|) (blue lines) and (ridi) 
(black lines) of level d by spins t and J, as a function of 
Vg, for the 1-orbital quantum dot circuit of Fig. Q The 
results are shown for the same parameters as in Fig. |5|and 
lead polarizations in the parallel configuration (c = P), with 
It-sdips = (panel a), gt^nhgijips ~ 0.06C/, (panel b) and 
gfJ-BSDiPS = 0.4(7 (panel c). For h^oips = 0, {ud]) and 
(ridi) remain very close, simply showing two steps correspond- 
ing to the two conductances peaks visible in Fig. |2Ja. In the 
case of a finite h^BiPS^ {''^di) rises more strongly than {udy) 
at the first conductance peak, revealing that current trans- 
port is due in majority to ] spins for this first peak. Then, 
(ndj) and {nd\) become closer when both £^d-\ and £_di are be- 
low the Fermi level. At the second conductance peak {jidy) 
rises more strongly than {nd-\) because current transport is 
now dominated by down spins. The asymmetry between the 
behaviors of spins t and i increases with h^oips (from left 
to right panels). 



FIG. 2: Panels a, c and e: Conductance in the parallel 
configuration (red full lines) and conductance G^^ in the an- 
tiparallel configuration (black dotted lines) as a function of 
the gate voltage Vg, for the circuit shown in Fig. Q with M 
a 1-orbital quantum dot. We have used Tl ~ 0.00517, Fp = 
0.07(7, Pt(fl) = 0.2, U/kpT = 30 and h^EiPs = 0. Panels b, 
d, and f: Magnetoresistance MR = (G^ - G'*^)/(G^ + G'*^) 
(pink curves) corresponding to the left conductance plots. 
The results are shown for g^J-ph'sDiPS = (panels a and b), 
gfJ-phg^upg — 0.06(7 (panels c and d) and gfJ-ph^DiPS ~ 
0.4(7 (panels e and f). 



corresponding to level d is split due to Coulomb inter- 
actions (see Eq. (gj, Figs. HJa, [2}c andEJe). At low 
temperatures T < Tk, Kondo effect is expected in the 
valley between the two resulting peaks. We have checked 
that the hypothesis T > Tk and hence the E.O.M. tech- 
nique are valid for the parameters of Fig. |2] (see Refs. 
lillia) . For hgDiPS = Oi we already note a strong quali- 
tative difference with the non-interacting case: although 
the two conductance peaks displayed by G^{Vg) are very 
similar, the MR variations corresponding to these two 
peaks have different shapes^. More precisely, for the low 
values of polarization considered here, MR(V^) is approx- 
imately mirror symmetric from one conductance peak to 
the other. Note that in Fig. 13 we have used specific 
parameters such that MR remains positive for any value 
of Vg when there is no SDIPS. Nevertheless, it is possible 
to have MR < for h^ujpg — h^DiPS ~ f*-"" instance 



by increasing Pl(R) {'c^oi shown). 

We now address the effect of a finite effective field 
hg^jpg. This field produces a shift of the conductance 
peaks from the P to the AP configurations. For instance, 
in Fig. 2-c and Fig. 2-e, plotted for h^pupg 7^ and 
^SDiPS ~ 0; the left [right] conductance peak is shifted 
to the right [left] from P to AP because it mainly comes 
from the transport of up [down] spins in the P case (this 
can be seen from the average occupations of the levels 
versus Vg in Fig. |3J). As a consequence, in Fig. [3 MR 
becomes negative for certain values of gate voltage. The 
effective field hgp,jpg thus enhances negative MR effects, 
ff hgDiPS i^ strong enough, it can even produce a giant 
MR effect with its sign tunable with Vg (Fig.HJf). More- 
over, because of the opposite shifts of the two consecutive 
conductance peaks for c = P with respect to those for 
c = AP, the positive-and-then-negative profile of MR 
corresponding to one conductance peak is generally fol- 
lowed by the negative-and-then-positive profile near the 
next conductance peak (approximately mirror symmet- 
ric). Note that a sign change h^p)ips ~^ ~^sdips '^^^^ 
not modify this behavior for the low values of polariza- 
tions considered here, because the spin with lower[higher] 
orbital energy will dominate in the left [right] peak of 

G^iVg). 

We now compare the results of this section with the ex- 
perimental data of Ref. Like many Coulomb blockade 
devices, the circuit studied in this experiment suffered 
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from low frequency T^g-noise, which can be attributed to 
charge fluctuators located in the vicinity of the device. A 
strong gate voltage offset jump occurred a.tVg = 4.331 V, 
and the data before and after this jump do not necessar- 
ily correspond to the filling of consecutive levels. There- 
fore, we will focus on the data taken for Vg > 4.331 V, 
shownd^ with black squares in Fig. IHl These data dis- 
play almost 2 regular MR(V^) oscillations, which cannot 
be understood with the 1-orbital model. Indeed, as ex- 
plained above, in this model, the two conductance peaks 
of C^Vg) are shifted in opposite directions by h^DiPS- 
As a consequence, the MR(V^) variations corresponding 
to these two peaks cannot be similar for parameters con- 
sistent with the experiment. We have shown here curves 
for hg£jpg ~ 0, but a finite hg^jpg would not modify 
this result. Using values of Pl(r) larger than in Fig. [21 
would not help either. 



For simplicity, we have considered in this section the 
one-orbital case. In reality, there are more than one or- 
bital levels on a quantum dot. As long as these levels are 
sufficiently well separated from each other (roughly, by 
an orbital energy difference larger than the Hund-rule ex- 
change energy), the two conductance peaks associated to 
a given level will occur consecutively in G'^{Vg) and will 
thus be described qualitatively like above^. In particu- 
lar, the two peaks will be shifted in opposite directions 
by h'^sDiPS peak to lower values of Vg and 

the second peak to higher values. Therefore, this limit 
should not allow to obtain two consecutive conductance 
peaks with analogue MR{Vg) patterns. On the contrary, 
if two (or more) levels are nearly degenerate, it is possible 
that the orbital levels of the quantum dot are not filled 
one by one while increasing Vg. Therefore, consecutive 
conductance peaks may exhibit a qualitatively different 
behavior compared with the one-orbital case. To examine 
this effect, we will consider in next section the extreme 
case of a quantum dot with a twofold orbital degeneracy. 
We will see that the discrepancy between the theory and 
the data can be resolved by using this model. 



Before concluding this section, we make a remark on 
another possible contribution to the spin splitting of the 
conductance peaks. Even though so far we have mainly 
considered the contribution from the SDIPS, in princi- 
ple, virtual particle exchange processes with the spin- 
polarized leads can also renormalize the energy levels 
through the ^Xd'a' ^^rms of Eq. ij^J^L^. Indeed, the 



^'dad'a' tcrms are not negligible in general. For exam- 
ple, in Fig. El-a, they globally shift the position of the 
conductance peaks in G^{Vg) by about 3.2% of U/ea. 
Nevertheless, for the low values of polarizations Pl{r) 
and the temperatures used here, the level spin-splitting 
produced by the S^'^^,^, terms is much weaker than this 
global shift and cannot compete with the finite values of 
hgujPS considered in this article. 



IV. QUANTUM DOT WITH A 
DOUBLY-DEGENERATE LEVEL 

In order to improve the understanding of Ref. 0, 
we now take into account the K-K orbital degeneracy 
commonly observed^'^i^^i^^i^^i'^^i^^ in SWNTs, by consid- 
ering a two-orbitals model i.e. hamiltonian ^ with d G 
{K, K } and S^K'a = £,Ka- Interestingly, SU(4) Kondo 
effect involving the orbital and spin degrees of freedom 
was observed in SWNTs with the K-K degeneracy^^'^°. 
This suggests that, in this system, the orbital quantum 
number is conserved during higher order tunnel events, 
probably because the electrons of the nanotube quantum 
dot are coupled to the nanotube section underneath the 
contacts, where they dwell for some time before mov- 
ing into the metal. For simplicity, we will also assume 
such a situation here and disregard high-order quantum 
processes which couple the K and K orbitals-^:. In or- 
der to calculate the conductance of the system from Eq. 
©, one needs to calculate the retarded Green's function 
Gs{u:) for s G {{K ]}, {K j}, {K' T}, {K' [}}. For this 
purpose, we again use the EOM technique. Since it is 
not possible to obtain a simple analytical expression for 
Gs{i-ij) in the two-orbitals case, we show below the system 
of equations of motion calculated by neglecting electronic 
correlations between the dot and the leads (T > Tk)- Us- 
ing s, si, S2 and S3 to denote four different dot states in 
the ensemble {{K ]}, {K j}, {K j}, {K i}}, we obtain 



Gs 



{hw - £,s 



{h+UiD^^+D^^+Dt^)} 

(12) 



Dt^ = (;jc.-6-c/-s°-sli)"'{^Ki>-Ss.iG. 

(13) 



,1,1 



(14) 



.3,5 



^S.Si ^S.,S2 ^S,S3/ 



(15) 



Due to interaction U, the Green's function 
Gs — Gs{uj) is coupled to other Green's func- 
tions Dfi'---'"^(w) = '•••''*3(w)e''^*dt with 

, - ~i9it){{n,,{t)..n,,it)cs{t),cl}) and 
nsi{t) — cl.{t)csi(t) for i e {1,3}. This means that 
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FIG. 4: Panels a, c and e: Conductance in the parallel 
configuration (red full lines) and conductance G^^ in the an- 
tiparallel configuration (black dotted lines), for the circuit of 
Fig. Q with M a two-orbitals quantum dot. We have used 
identical tunnel rates to the two orbitals, i.e. Fz, = 0.0043?7, 



Tr = 0.0725f/, and Pl(r) = 0.4. We have also used 
U/knT = 30 and hj^jps = 0. Panels b, d and f: Mag- 
netoresistance MR (pink full lines) corresponding to the left 
conductance plots. The results are shown for g^-sh^DiPS = 
(panels a and b), gfJ-sh^j^jpg = 0.05f/ (panels c and d) and 
gUBhgj^ijpg — 0.3U (panels e and f). 




FIG. 5: Average occupations (n|) — (uki) = {?^i^'t) (blue 
lines) and (n;) — {tiki) = (nK'i) (black lines) of levels K 
and K' by spins | and J, as a function of Vg, for a two- 
orbitals quantum dot circuit with the same parameters as 
in Fig. |1] The results are shown for lead polarizations in 
the parallel configuration (c = P), with h^oips = (panel 
a), g/J-Bhgpijpg — 0.05U (panel b), and g^,Bh^DiPS = 0.317 
(panel c). For h^oips ~ 0, {n-^) and {n^) remains very close, 
showing four steps corresponding to the four conductances 
peaks visible in Fig. In the case of a finite h^oips^ {'T-t) 
rises more strongly than {n{) for the two first conductance 
peaks, which shows that current transport is due in majority 
to up spins for these two peaks. The opposite situation occurs 
for the two last conductance peaks. The asymmetry between 
the behaviors of spins | and J, increases with h^oips (from 
left to right panels). 



Eq.®, and terms defined by 
/ 



^da,d'i 
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^d' 



^ -idu - id' a- + ika' - nil + ^0^ 
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fJ-tiika') 


^d' 


2 


— idcr — id' a' + Cfecj 




n + 2)U + iQ+ 



(16) 



the dynamics of electrons in state s is modified by the 
presence of other electrons on the dot [In the one orbital 
case, Gdcr was coupled to only, which lead to simple 
expression ©]. The term E° = -i{T^ + Ff )/2 is the 
tunneling self energy for a non-interacting quantum dot, 
already introduced in previous section. The equations of 
motion also involve terms S^'"^,^, , which are defined by 



These terms take into account the tunneling of elec- 
trons between the leads and a dot state d'a' differ- 
ent from da. The average level occupations occur- 
ing in Eqs. lfT^ - ((T5|) are given by (rtgj , .., n^^) — 
- d^f{huj)lm[Dn'-''''^-'{uj)]. The Green's func- 
tions and the level occupations can be calculated numer- 
ically from the above equations. Up to now, the E.O.M. 
technique for multilevel systems had been implemented 
only by neglecting S*^'"_^,^, and Xda,d'a' termsS?. Like in 
the one-orbital case, these terms are not negligible in the 
context of our study. 

Figure 0] shows the conductance (panels a, c and e) 
and MR curves (panels b, d and f) calculated from Eqs. 
Q and H12I15|I . for different values of hf^gjjips- For sim- 
plicity, we have assumed that the coupling to the leads 



is identical for the two orbitals, i.e. T^^^j^j = r^'|j|j = 
ri(l±P/) for I e {L, R} . We have again used parameters 
consistent with Ref. [13 i.e. U /ksT = 30, relatively low 
polarizations \ Pl{r) \ = 0.4 and values of ri(/j) leading to 
the proper width and height for the conductance peaks. 
We have checked that these parameters are compatible 
with the hypothesis T > Tk, with Tk the Kondo tem- 
perature associated to the SU(4) Kondo effect expected 
in this system^. In most cases, the curves G^(V^) show 
4 resonances, the first two associated with a single occu- 
pation of K and K', and the other two to double occupa- 
tion (see e.g. Fig. ^a.). For hg^jpg = hg£jpg = and 
the parameters used here, MR remains positive for any 
value of Vg (Fig. 01-b). Like in the 1-orbital case, a finite 
^SDiPS makes easier negative MR effects and can even 
lead to a giant MR effect with a sign tunable with Vg 
(Figs 4-d and 4-f). Importantly, the effect of h'^sDiPS 

de- 
pends on the occupation of the dot. For instance, in Fig. 
21-e plotted for gfJ'shgjjjpg larger than the linewidth of 
the conductance peaks, the first two conductance peaks 
of (peaks 1 and 2) are strongly shifted to the left by 
^SDiPS because they are due in majority to up spins, as 
can be seen from the average occupation of the levels in 
Fig. 13-c. This allows to get a MR pattern approximately 
similar for these two peaks, i.e. a transition from posi- 
tive to negative values of MR (Fig. 0}f ). On the contrary, 
peak 4 corresponds to a transition from negative to pos- 
itive values of MR because the associated conductance 
peak is due in majority to down spins. In Fig. ^J-e, the 
shape of the MR(V^) pattern associated to peak 3 is more 
particular (positive/negative/positive) because, for the 
values of parameters considered here. Coulomb blockade 
does not entirely suppress the up spins contribution in 
peak 3, which is therefore spin-split-^'' . Remarkably, this 
allows to obtain, at the left of Fig. 0J-f, three positive MR 
maxima which differ in amplitude but have rather simi- 
lar shapes. In the case of gi^'sh^DiPS finite but smaller 
than the linewidth of the conductance peaks (Fig. El-c), 
the amplitude of the MR signal is much smaller than in 
the previous case but its shape remains comparable. 

We now reconsider the experimental data of Ref. IT^ 
Even the two-orbital model cannot not provide a rea- 
sonable fit to the data if we assume h^jjjpg = and 
^SDiPS ~ ^- contrast, the two-orbital model ex- 
hibits a good agreement with the experimental data for 
OmU, h^'oips = 0, Tl/U = 0.0043, Tr/U 
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FIG. 6: Comparison between the data of Ref. ^3 (squares) 
and the two-orbitals theory. We show the conductance in 
the parallel configuration (top panel) and the corresponding 
magnetoresistance MR' = - G^^ ) / G"^^ (bottom panel). 
The theory is shown for parameters consistent with the ex- 
periment i.e. ?7 = 5 meV, U/kBT = 30 and a = 0.0986. 
We also use relatively low values of polarization Pl(r) ~ 0.4 
because usual ferromagnetic contact materials are not fully 
polarized^. Assuming identical tunnel couplings for the two 
orbitals, the values of tunnels rates Tl ~ 0.0043?7 and Tr = 
0.0725(7 are imposed by the width and height of the conduc- 
tance peaks. Then, hg\^j^pg are the only free fitting param- 
eters which remain for interpreting the MR curve. We have 
assumed g^sh^DiPS ~ 0.05(7 for all the theoretical curves 
shown in this Figure. We have plot the MR curves of bottom 
panel for h^pjps ~ (pink curve, corresponding to Fig. |ljd), 
h^BjPS ~ —0.01(7 (green dashed curve) and h^pjps = 0.01(7 
(blue dot-dashed curve). Note that in this Figure, we show 



MR' instead of MR = (G^ - G'*^)/(G^ + G^^ 



in order 



"-SDIPS 



to be consistent with Ref. II 7l . A strong gate voltage offset 
jump occurred at Vg = 4.331 V, therefore, we show the data 
at the left/right of this jump with grey/black symbols. 



estimated /isd/ps ^^e non- interacting theory us- 
ing realistic parameters i.e leads with a Fermi energy 
10 eV and a density of states polarized by 40%, and 
a nanotube with Fermi wavevector 8.5 lO^rn^^, Fermi 

M 



Ref. [13, and density of states N^' = 2L/nhvP. We have 
modeled the interfaces between the nanotube and the 
leads with Dirac potential barriers^, with a height which 
is spin-polarized by 40% and an average value which cor- 
responds tc^^ Yp^p^ = Tl(r)/2ttNp'-^^ ~ 60 /zeV (For 
comparison the fitting parameters used in Fig. corre- 
spond to Fi = 21 ^eV and Tn = 362 ^eV). We obtain 
^SDiPS ^ 1-3 T, which is consistent with the above anal- 
ysis. 

In the above discussion, we have assumed h^oips ~ 
for simplicity. The height of the conductance peaks in 



0.0725, \Pl(R)\ = 0.4, and parameters U = b meV, velocity = 8 lO^m.s \ length L = 500 nm like in 
U/kBT — 30, and a = 0.0986 given by the experiment 
(see Fig. 1^1 red and pink full curves). 

We now discuss the value of hgpjpg — 0.05(7 found 
for the above fit. This corresponds to a magnetic field 
of about 2 T, which is too strong to be attributed to 
stray fields from the ferromagnetic electrodes (see e.g. 
Ref. lU). This is in favor of generalizing the SDIPS con- 
cept to SWNTs quantum dots circuits, i.e. considering 
that the energy levels of the dot are spin-split because 
the confinement potential created by the ferromagnetic 
electrodes is spin-dependent. For comparison, we have 
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the data imposes to use a strong assymetry Tji/rL ~ 
17 between the left and right tunnel rates. Thus, the 
two tunnel barriers are not symmetric, and there is no 
fundamental reason to assume hg£jpg = 0. Figure |^ 
shows examples of MR curves plotted for a finite h^DiPS- 
Using hg^jpg = 0.01J7 (green dashed curve) enhances 
the fit of the MR at peak 2 whereas h^uips ~ — O.Olt/ 
(blue dot-dashed curve) enhances the fit of the MR at 

peak 3. Interestingly, with the non-interacting model, 

p x P \ 

assuming the most simple situation in which (/Sj — 
has the same sign for the two leads, one finds j/i^lj/pgl < 
I ^SDiPS I ' which is in agreement with the values used 
here. The fact that the best fit for the MR patterns at 
peaks 2 and 3 correspond to different values of hg£jpg 
might be due to a gate dependence of the SDIPS. This is 
indeed possible since the potential profile of the interfaces 
between the wire and the leads can vary with Vg. 

We now comment briefly on the data taken for Vg < 
4.331 V. It is not sure that the data shown before and 
after Vg — 4.331 V correspond to the filhng of consecutive 
levels because of the gate voltage jump which occured at 
this value of Vg. Nevertheless, the shape of the MR curve 
corresponding to Vg < 4.331 V is rather consistent with 
the theory shown in Fig. 4. This suggests that these 
data really correspond to peak 1. At this stage, it is 
important to point out that other orbital levels not taken 
into account in our calculation should slightly modify the 
conductance peaks 1 and 4. The discrepancy between the 
theory and the data for Vg < 4.331 V could be explained 
by the effect of the other orbitals. For the data at Vg > 
4.331 V, our fit is more quantitative since we have used 
peaks 2 and 3 of the theory. 

In principle, the modelisation of the orbital levels in 
SWNTs can be refined by taking into account an ex- 
change energy J which favors spin alignment, an ex- 
cess Coulomb energy 6U related to the double occu- 
pation of the same orbital, and a subband mismatch 
<^ = ^K'a — £,Ka 7^ (see Ref. ISTj) . In practice, SU 
is rather small but J and S can be of the same order 
as U. Two different regimes of parameters can occur in 
practice, li 6 > J + SU + Igf^phgoipsl' electrons 
with opposite spins will fill consecutive ly the same en- 
ergy level while Vg increases (see Refs. 20 30), and the 
behavior of the device should thus be analogue to the 
non-degenerate multi-orbital case evoked at the end of 



section UTTI Nevertheless, ii S < J + 5U + Isms'isd/psI' 
peaks 1 and 2 [3 and 4] will correspond in majority to the 
same spin direction, as observed experimentally bv l23l2^ 
In this case, the effect of h'^sDiPS should be qualitatively 
the same as described in the present section. We expect 
that the weights of K and K' in peaks 1 and 2 differ due 
to (5 7^ 0, but this should not change the way in which 
^SDiPS shifts the conductance peaks from P to AP. 

In future experiments, it would be interesting to obtain 
continuous data on a larger Vg-range, in order to check 
that the shape of the MR{Vg) pattern depends on the 
occupation of the dot. This would also allow to study the 
gate voltage dependence of h'^sDiPS- would also be in- 
teresting to engineer contacts with ferromagnetic insula- 
tors or highly polarized ferromagnets in order to observe 
the SDIPS-induced giant MR effect. Note that although 
we have considered here the limit /c^T < F^^ -I- F^^., 
a strong enough SDIPS should also affect the behav- 
ior of the quantum dot in the sequential tunneling limit 
/cfiT » F^^ + F^^, through an analogous mechanism. 

V. CONCLUSION 

Using an Anderson model, we have studied the behav- 
ior of a quantum dot connected to ferromagnetic leads 
through spin-active interfaces. The spin activity of the 
interfaces makes easier negative magnetoresistance (MR) 
effects and can even lead to a giant MR with a sign oscil- 
lating with the gate voltage of the dot. Due to Coulomb 
blockade, the MR versus gate voltage pattern cannot be 
identical for all conductance peaks. It is nevertheless pos- 
sible to account for the MR data measured by Ref. 11 in 
single-wall carbon nanotubes by taking into account the 
K — K' orbital degeneracy commonly observed in those 
systems. 
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